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1. INTRODUCTION 


The major objectives of our lunar studies have been to determine the origin 
and history of the field that produced the remanent magnetism on the moon, 
and to obtain Information on the history and present state of the lunar 
crust. In previous years we have concentrated on the production of maps of 
the lunar magnetic field and Interpretation of the data represented by the 
maps. During the past year, we studied various approaches to the determina- 
tion of mathmatical representations of planetary magnetic fields fr«n 
satellite surveys of the fields. 

The main emphasis of our research effort during the past year has been on 
determining methods for the optimal processing of satellite-derived marnetlc 
anomaly data. The motivation for undertaking this effort was not simply the 
expectation of being able to make modest Improvements In our lunar maps, but 
rather the expectation of being able to produce vsstly superior maps and 
models of the lunar magnetic field with respect to accuracy and spatial 
coverage, and of being able to describe this field statistically for compari- 
son with the gecT>agnet1c field. 

The longer-term objective Is knowledge of the characteristics of the Induced 
and remanent magnetism of bodies in the solar systas as well as the plane- 
tary dynamo magnetic fields of these bodies. 

(^estlons of a fundamental nature exist regarding the Moon's crustal magnetic 
field, ftost importantly, the origin of the field that produced the crustal 
magnetization temains unknown. Specification of its origin would have impli- 
cations for the early history of the solar system as »«11 as that of Moon. 
Although near-surface magnetic field measurements are available only for 
Earth and Moon, it seens likely that crustal magnetization is characteristic 
of all the silicate planetary bodies in the solar system. Thus the develop-- 
ment of methods for the quantitative description of crustal magnetic fields 
is a task of significance for future planetary exploration. The modelling 
of the lunar field with equivalent point sources, as suggested herein, is a 
way to succinctly describe the crustal lunar field. The magnetic field at 
various altitudes can be readily computed from such a model and represented 
by contour maps. Thus this model is an important first step toward the ulti- 
mate goal of understanding the processes that have produced the crustal 
lunar field. 

With respect to future lunar and planetary missions, the results of our roost 
recent studies can be used to establish requirements for amplitude and tem- 
poral resolutions of the field measurements, establish requirements for the 
altitude range of the orbits, and identify regions of special interest. 

2. LIMITATIONS OF PRESENT LUNAR KAGNETIC MAPS 

Correcting for altitude differences among the measurement positions when con- 
structing a magnetic anomaly map from satellite-derived data is normally a 
nontrivial problem (Bhattacharyya, 1977; Mayhew, 1979; Mayhew et al., 1980). 
The slow lunar rotation period allows a simpler solution for a subset of the 
Apollo data. If only data collected during a single lunation over a par- 
ticular stenographic area are used, the magnetic anomaly map can be defined 
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on the curved surface of observation of the subsatellite. This approach 
Mas followed by Hood et a1. (1981). However, It has a nwnber of serious 
limitations. 

a. Limited selenographic extent. Although data were collected between 
t30** latitude, only about fifteen percent of this region has been 
mapped. 

b. Small fraction of available lunar data used. Even for the limited 
region mapped, only a small fraction of the available data has been 
used. 

c. Maps not for spherical surface. The field was mapped on the sur- 
face defined by the spacecraft orbits. This surface Is "wavey" 
and of varying altitude, and not convenient for upward or downward 
continuation. 

d. Filtering not optimized. 

e. Error analysis not inherent in the procedure. 

f. Data not in convenient form for statistical studies. 

The limitations noted are a direct consequence of the n^thod of data pro- 
cessing chosen. It is nearly impossible to overstate the severity of these 
limitations; in comparison to what is potentially achievable with respect 
to accuracy and spatial coverage of lunar magnetic maps and models based on 
the Apollo data sets, the maps published thus far represent only a crude 
and primitive beginning. 

Because we recognized the serious limitations of these lunar magnetic maps, 
v« initiated a study during 1980 with the objective of determining methods 
for the optimal processing of satellite-derived magnetic anomaly data. The 
results of this study represent a considerable advance beyond the methods 
for planetary magnetic field mapping and modelling that were available at 
the time our production of lunar maps was initiated. The methods available 
at that time were essentially the ones used on the first geomagnetic ancwialy 
map, based on POGO satellite data, produced by Regan et al. (1975). 

3. OPTIMAL PROCESSING OF SATELLITE-DERIVED MAGNETIC ANOMALY DATA 

In the course of a design study concerned with satellite mapping of the crus- 
tal magnetic field of the Earth, we undertook to investigate methods by which 
to optimize the production of a mathematical representation of the field fr«n 
a dataset acquired with orbiting spacecraft. During a previous reporting 
period, we initiated an examination of the similarities and differences of 
the problems of mapping the crustal fields of the Earth and Moon. 

Three papers pertinent to both problems, and thus to the general problem of 
optimal mapping of planetary magnetic fields from orbiting spacecraft, were 
published during that previous reporting period. In one, McLeod and Coleman 
(1980a) discussed some of the limitations of geomagnetic field mapping in 
terms of the power spectra of instrumentation "noise'' and of "noise" due to 
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time varying external current systems relative to the power spectra of the 
"signal” due to spacecraft motion through the spatially varying geomagnetic 
field. It was shown how these spectra are related to the achievable resolu- 
tion for a given spacecraft altitude, and how these spectra influence the 
optimal data processing procedure. The second paper, McLeod and Coleimn 
(1980b), is concerned with the statistical description of the geomagnetic 
field by means of spherical harmonic power spectra, with numerical values 
of these spectra for the Earth, and with statistical field models that would 
produce the observed spectra. This paper also contains a discussion of the 
relation between spherical harmonic power spectra and great circle Fourier 
power spectra for a circular orbit. In the third paper, McLeod (1980) de- 
rives orthogonality relations used in McLeod and Coleman (1980b). »These 
papers are relevant not only to geomagnetic field mapping, but also to map- 
ping the magnetic fields of the moc;; ind other planets. 

During that previous reporting period, we presented a paper based on these 
three publications at the NASA Geodynamics Program review at GSFC, January 
26-29, 1981 (McLeod and Coleman, 1981). 

During the current reporting period, we continued our studies of planetary mag- 
netic field nwdelling and mapping. We anticipate that the results of these 
efforts will find application not only in processing satellite-derived mag- 
netic anomaly data for the Moon, Earth, and other planets, but also in pro- 
cessing aeromagnetic data (for Earth) of much shorter wavelength. 

A paper detailing the results of these studies has been completed and ac- 
cepted for publication (McLeod, 1982a). The paper contains an analysis of 
the data processing method used by Mayhew et al. (1980) on POGO satellite 
data over Australia. It is shown that the data processing method used by 
Mayhew et al. is approximately equivalent to filtering the data to retain 
only a specific region of the two-dimensional frequency plane. The filter 
is reasonably close to optimal with respect to accuracy, though some pos- 
sible improvements are suggested. The analysis makes use of the concept of 
two-dimensional spatial power spectra and the ideas of Wiener (1949) as dis- 
cussed by Davenport and Root (1958). Errors remaining after the data pro- 
cessing are discussed, and the analysis indicates how to estimate these 
errors. The techniques used can be carried over to the problem of lunar 
magnetic field modelling and mapping. The studies indicate that the basic 
principles used by Mayhew et al. for crustal geomagnetic field modelling and 
mapping appear to be most efficient for modelling and mapping the lunar 
field; however, a number of changes in the detailed methods are necessary 
and desirable. 

A direct harmonic analysis of the data would be impractical for processing 
the data obtained from a low altitude satellite such as Apollo or Magsat. 
Harmonic analysis, using methods discussed by McLeod (1982a), can be a use- 
ful technique for statistical studies of the lunar magnetic field. 

During this reporting period, Dr. McLeod attended two scientific meetings 
devoted to geodynamics and geomagnetism and made presentations related to 
the optimal processing of satellite-derived magnetic anomaly data at each 
meeting. Abstracts of the presentations (McLeod, 1982c, d) were published 
in the proceedings of these meetings. An abstract (McLeod, 1982e) for a 
presentation scheduled to be made at the spring AGU meeting in Philadelphia 
will be published in the May 4, 1982 issue of EOS. 
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4 . STUDIES OF CRUSTAL AND CORE GEOMAGNETISM 

Two other papers completed during this reporting period (McLeod, 1981, 1982b) 
deal with the general subject of geomagnetism as related to both the core and 
the crust. Besides the Imnedlate scientific value of these studies, the 
mathematical methods employed are relevant to the study of other planetary 
fields. Including the lunar field. 

The first paper (McLeod, 1981) Is concerned with geomagnetic variations hav- 
ing (apparent) periods from 13 to 30 years as reported by Alldredge (1977), 
who has argued that the origin of these signals must be found In the core 
of the Earth rather than outside the Earth. It Is shown In this paper that 
a portion of these geomagnetic variations (perhaps most of the variations) 
might well be due to geomagnetic signals of much longer period, originating 
In the core, that appear to have periods of 13 to 30 years because of an arti- 
fact of the data processing. Much of the remainder of these variations could 
well be of external origin. A superior method for processing these data is 
suggested in the paper. 

Also In this first paper, the background spectra of the geomagnetic field 
(I.e., with "lines" In the spectra removed) reported by Courtlllot and Le 
Mou^l (’976) for the period range 10 years to 2 ninths are discussed. These 
spectra are approximately of the (1/f') type, as are the geomagnetic spectra 
for the vastly different period range 4 hours to 5 minutes reported by Camp- 
bell (1976) and discussed by McLeod and Coleman (1980b). It is noted that 
If Campbell's spectra are extrapolated to the period range investigated by 
Courtlllot and Le Mou^l, an extrapolation over nearly 3 decades, the two 
spectra agree numerically (approximately). It Is further noted that this 
(l/f^) spectrum Is the spectrum corresponding to the well-known random walk, 
or "drunkard's walk." 

The second paper (McLeod, 1982b) is concerned with the crustal geomagnetic 
field and the statistical description of this field by n«ans of spatial 
power spectra. Two-dimensional Fourier spatial power spectra of equivalent 
magnetization values for a region that includes a large portion of the 
western United States are presented. The magnetization values were deter- 
mined by Inversion of POGO satellite data, assuming a 40 km thick magnetic 
crust, and were located on an 11x10 array, 300 krr grid spacing. The spectra 
appear to be in good agreement with values of the crustal geomagnetic field 
spatial power spectra given by McLeod and Coleman (1980b) and with the crus- 
tal field model given by Serson and Hannaford (1957). The spectra show evi- 
dence of noise at low frequencies in the direction along the satellite 
orbital track (N-S), indicating that for this particular data set, addi- 
tional filtering would probably be desirable. These findings Illustrate the 
value of two-dimensional spatial power spectra both for describing the geo- 
magnetic field statistically and for diagnosing possible noise sources. 

5. LUNAR REMNANT MAGNETIZATION 

During a previous reporting period, the analysis of Apollo 15 and 16 subsatel- 
lite magnetometer data was extended as one element of a study of the nature and 
origin of lunar crustal magnetization. This work was supported with funds 
from NASA grant NS6 7020 LS, the University of Arizona, C. P. Sonett, Prin- 
cipal Investigator, as well as funds from this grant, and performed In 
collaboration with L. L. Hood at Arizona (Hood et al., 1981). 
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Correlations of mapped anomalies (Hood et al.« 1981a) with surface geologic 
units continue to Indicate that the n»st strongly magnetized materials may 
not be deep seated In the crust as expected If their magnetization was ac> 
quired by slow cooling In the presence of a steady magnetic field. The least 
ambiguous correlations are found for relatively surf 1 dal units such as the 
Fra Mauro Formation and the ftelner Ganma swirls (Hood» 1980b). A second 
example of a magnetic anomaly maximum detected wUh the subsatellite mag- 
netometers over an area dominated by swirl markings was found near crater 
Ellerman. Although the radial field amplitude Is less than 2 gammas (1 gamma 
* 10"5 G) at 87 km above the surface, an approximate scaling for altitude 
differences Indicates that the source dipole moment Is comparable to or 
larger than that of other anomalies detected with both subsatellites Includ- 
ing the Reiner Gamma anomaly and the Van de Graaff-Aitken anomalies. This 
concentration of swirls is located within 5® of the Crisium antipode raising 
to at least four the proposed number of young large Impact basins that are 
nearly antipodal to concentrations of swirls and strong magnetic ancxnalles 
detected by both magnetometers and charged-particle Instruments on board 
the Apollo subsatellites. 

The Frau Mauro Formation Is Imbrium basin ejecta. There is Independent evi- 
dence that areas peripheral to other large impact basins are also magnetized. 

If the swirls are produced by either meteroid or cometary Impacts as has been 
suggested, then It remains possible that all observed anomalies are associ- 
ated with Impact-generated surface materials. 

A method for distinguishing between local surface magnetizing fields such as 
may have been produced during impacts and a large-scale or global magnetiz- 
ing field (Srnka, 1977; Gold and Soter, 1976) and a large-scale or global 
magnetizing field Involves estimating the directions of the bulk magnetization 
of relatively strongly magnetized and isolated magnetic anomaly sources using 
the subsatellite vector magnetometer data. In a previous study (Hood et al., 
1978) that was limited primarily to the Van de Graaf f-AI tken region on the 
south-central far side. It was found (a) nearby sources are generally magne- 
tized in very different directions; (b) the number of north-south oriented 
vectors is depleted relative to the number of east-west and radially oriented 
vectors; and (c) the inferred magnetization directions are otherwise nearly 
random. Assuming that these properties were repeated elsewhere on the moon, 
several interpretations were suggested. First, the negligible lunar mag- 
netic dipole moment (Russell et al., 1974) was explained as due to the In- 
ferred lack of directional coherence for the magnetization rather than as a 
consequence of global crustal magnetization along dipolar field lines as ear- 
lier suggested (Runcorn, 1975). Second, the north-south depletion of the 
magnetization Intensity was considered to represent an observational constraint 
that did not favor purely local surface origins of the magnetizing field. An 
exception was meteorltic or cometary Impact processes that might tend to am- 
plify the weak solar wind magnetic field which, on the average. Is parallel 
to the lunar equatorial plane. 

During the previous reporting period, additional maps of the vector components 
of the crustal magnetic field became available with low-altitude coveraoe ex- 
ceeding that of older maps by a factor of about ten (Hood et al., 1981a;. 
However, even these maps are of very limited selenographic extent, covering 
only about fifteen percent of the region for which data are available, as 


Indicated in a previous section. Thus further Improvements In selenographic 
mapping and modelling would be very valuable. 
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Bulk directions of magnetization were estimated for sources of the strongest 
and most Isolated magnetic anomalies apparently associated with exposed rem- 
nants of the Fra Mauro Formation, and a consolidation of swirls near crater 
Ellerman. A total of twenty-five directions of magnetization excluding the 
Van de Graaff-Altken region have thus far been estimated for sources that 
are widely distributed in longitude but that are confined latitudinally by 
the near-equatorial orbits of the Apollo subsatellites. The results again 
show that adjacent sources are. In general, magnetized In very different 
directions. This has already been documented for the Important special case 
of the Fra Mauro anomalies (Hood, 1980b). A further result Is that the pole 
positions of a magnetizing dipole assumed to be centered in the moon cor- 
respond to these twenty-five vectors are not clustered In any discernible 
manner. This finding contrasts with that of Runcorn (1979) who found a sig- 
nificant clustering of pole positions corresponding to the bulk directions of 
magnetization inferred for the Van de Graaf-Aitken region. These results 
were presented to the Twelfth Lunar and Planetary Science Conference (Hood 
et al.. 1981b). 

In addition to the apparent absence of directional evidence for a global 
lunar magnetizing field, there is now the constraint imposed by the relatively 
small maximum size of the lunar metallic core radius as estimated from a vari- 
ety of geophysical measurements. In particular, the most recent upper limit 
from magnetic sounding data is about 360 km (Hood et al., 1981b). To produce 
a magnetic field at the lunar surface with an amplitude of order 1 G as re- 
quired by returned sample studies (14), the field at the surface of a 360 km 
core must exceed 100 G. For comparison, the corresponding field at the sur- 
face of the terrestrial core Is about 5 G. Final evidence from orbital data 
for magnetic fields generated by impact processes would cone in the form of 
a clear association of a mapped anomaly with a young (« 3 b.y. old) surface 
geological unit. Reiner Gamma is an obvious candidate but photogeologic 
considerations and the tendency for other swirls to be concentrated near the 
antipodes of large, young impact basins require a cautious Interpretation of 
these features (Hood, 1980a). 

6. INDUCED LUNAR MAGNETIZATION 

During a previous reporting period, Apollo 15 and 16 subsatellite fluxgate 
magnetometer data were analyzed for all Intervals In which the moon was in 
the lobes of the geo.nagnetic tail to obtain an improved estimate of the 
average magnitude of the induced dipole moment of the moon. Non tail lobe 
Intervals were excluded by an examination of particle data recorded simul- 
taneously on the subsatellite and on the moon and by a comparison of the 
magnetic field strength observed with that expected for the solar wind condi- 
tions prevailing during the measurements. The surface magnetometer data were 
used to minimize effects due to temporal changes In the external field. The 
resulting set of estimates, while smaller than previous sets, has a smaller 
spread and yields an Induced magnetic moment of -4.23 x 10^^ Gauss-cm^ per 
Gauss of applied field, cor res ponding to a G-factor of 0.08 ± .011. 
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These results, reported to the Twelfth Lunar and Planetary Science Conference, 
do not place strong constraints on the conductivity of the lunar core. We 
would detect the observed effects as long as the core conductivity was greater 
than about 10 mho/tn. If the outer, cooler layers of the moon, that are at 
temperatures below the effective Curie point, contain little or no free iron, 
thetj these measurements are consistent with the presence of a conducting core 
whose radius is si ghtly larger than 400 km. If these outer layers of the 
moon contain significant amounts of free iron, and hence exhibit the para- 
magnetism expected in such a situation, the core si.Te could be ever greater. 

7. RECOMMENDATION FOR FURTHER RESEARCH 

The study of the optimal processing methods for satellite-derived magnetic 
anomaly data was motivated by the existence of serious limitations in the 
existing maps of the lunar field as published by Hood et al. (1981). These 
limitations are a result of the method of data processing used, so that while 
data were collected between i30® latitude, only about fifteen percent of this 
area has been mapped. The maps are not for a spherical surface but for an un- 
dulating surface defined by the spaceciaft orbits for a single lunation. It 
is believed that use of the methods described by McLeod (1982a) could result 
in an improvement in these maps by about an order of magnitude in both ac- 
curacy and spatial coverage. The resulting maps would be for a spherical 
surface, convenient for use in statistical studies, and would constitute a 
major improvement over existing maps and models of the lunar field. 

The general plan of the work should be as follows: 

1. Survey the data coverage of the UCLA Apollo magnetometer measurements. 

2. Filter the data by (partial) orbits to remove the lowest frequencies. 
Only data from the lunar wake should be used. 

3. Compute sets of normal equations arising from the least squares 
problem associated with fitting the measured field with the field 
due to equivalent point monopole sources, as discussed by McLeod 
(1982a). Northbound and southbound orbits should be treated sepa- 
rately. 

4. Solve the normal equations for various regions and combine results 
with attention to overlapping. 

5. Compute the resulting distribution of monopole sources on a finer 
grid with spatial filtering for input to contouring program. 

6. Produce various two-dimensional maps from these results. 

7. Subtract the computed values from selected samples of measured orbi- 
tal data to determine measurement error; compute power spectra of 
the error and use this computation to estimate errors of the previ- 
ously computed maps. 

8. Compute the spherical harmonic spectrum of the resulting source dis- 
tributions. 
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Some experimentation would be necessary to determine the optimal point source 

grid spacing, and the optimal filtering for steps 2 and 5. 
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